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We give a theoretical framework to obtain a low-energy effective theory of quantum chromodynam- 
ics (QCD) towards a first-principle derivation of confinement/deconfinement and chiral-symmetry 
breaking/restoration crossover transitions. In fact, we demonstrate that an effective theory obtained 
using simple but non-trivial approximations within this framework enables us to treat both tran- 
sitions simultaneously on equal footing. A resulting effective theory is regarded as a modified and 
improved version of nonlocal Polyakov-loop extended Nambu-Jona-Lasinio (nonlocal PNJL) mod- 
els proposed recently by Hell, Rossner, Cristoforetti and Weise, and Sasaki, Friman and Redlich, 
extending the original (local) PNJL model by Fukushima and others. A novel feature is that the 
nonlocal NJL coupling depends explicitly on the temperature and Polyakov loop, which affects the 
entanglement between confinement and chiral symmetry breaking, together with the cross term in- 
troduced through the covariant derivative in the quark sector considered in the conventional PNJL 
model. The chiral symmetry breaking/restoration transition is controlled by the nonlocal NJL in- 
teraction, while the confinement/deconfinement transition in the pure gluon sector is specified by 
the nonperturbative effective potential for the Polyakov loop obtained recently by Braun, Gies, 
Marhauser and Pawlowski. The basic ingredients are a reformulation of QCD based on new vari- 
ables and the flow equation of the Wetterich type in the Wilsonian renormalization group. This 
framework can be applied to investigate the QCD phase diagram at finite temperature and density. 

PACS numbers: 12.38.Aw, 05.10.Cc, ll.lO.Wx 



I. INTRODUCTION 



The relation betvifeen confinement and chiral symmetry 
breaking is one of the long-standing puzzles in theoreti- 
cal physics. Recently, strong interest on this issue revived 
in extreme environments especially at high temperatures 
and baryon densities, stimulated by the heavy-ion pro- 
grams at GSI, CERN SPS, RHIC and LHC, see e.g., 
[l|, |3] for a review. Quantum chromodynamics (QCD) 
for strong interactions is a fundamental theory for solv- 
ing this problem. 

In pure Yang-Mills theory, i.e., in the limit of infinitely 
heavy quark mass mq — )■ cx) of QCD, the Polyakov loop 
average {L), i.e., the vacuum expectation value of the 
Polyakov loop operator L, can be used as a criterion for 
quark confinement [3| . The Polyakov loop operator L is a 
gauge invariant operator charged under the center group 
Z{Nc) of the color gauge group SU{Nc). The Polyakov 
loop average (L) vanishes (L) — and quarks are con- 
fined at low temperatures T < Td where the global center 
symmetry Z{Nc) is intact, while it is nonzero (L) ^ and 
quarks are deconfined at high temperature T > Td where 
the global center symmetry Z{Nc) is spontaneously bro- 
ken. Thus, we can define Td as a critical temperature for 
confinement/deconfinement phase transition. 

When dynamical quarks in the fundamental represen- 
tation of the gauge group are added to the Yang-Mills 
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theory, the center symmetry is no longer exact. On the 
other hand, QCD with massless quarks rriq — >■ exhibits 
chiral symmetry SU{N-f)L x SU{Nf)R. The chiral con- 
densate (ipip), i.e., the vacuum expectation value of a 
gauge-invariant composite operator ip'ip, is used as an or- 
der parameter for chiral symmetry breaking. The chiral 
condensate (V'V') is nonzero (iptp) ?^ at low tempera- 
tures T < T^ where the chiral symmetry is spontaneously 
broken, while it vanishes (il'tp) = at high temperature 
T > Ty- where the chiral symmetry is restored. Thus, 
we can define T^ as a critical temperature for chiral- 
symmetry breaking/restoration phase transition. 

For realistic quark mass (with finite and nonzero ruq: 
< ruq < oo), there are no exact symmetries directly 
related to the phase transitions, since both the center 
and chiral symmetries are explicitly broken, and (L) and 
[ipip) are approximate order parameters. In this case, 
there is no critical temperature Tc in strict sense, and 
the transition can be a crossover transition for which the 
pseudo critical temperature T* is defined such that the 
susceptibility takes the maximal value at T = T*. If 
quarks are in the fundamental representation, deconfine- 
ment (a rise in the Polyakov loop average) happens at 
the temperature where the chiral symmetry is restored 
(chiral condensate decreases rapidly). The chiral and de- 
confinement transitions seems to coincide, T^ = T* c^ Tc 
d, 0], although the property of the phase transition, 
e.g., the critical temperature and the order of the tran- 
sition depend on numbers of color Nc and fiavor Nf. 
Whereas, for quarks in the adjoint representation, decon- 
finement and chiral-symmetry restoration do not happen 
at the same temperature, rather T* ^ T^ 0. Although 



there exist theoretical considerations on the interplay be- 
tween chiral symmetry breaking and confinement at zero 
baryon density [3] , the underlying reasons for the coinci- 
dence are still unknown and uncertain at nonzero baryon 
density [M3. 

The hadronic properties, especially, chiral dynamics 
at low energy have been successfully described by chi- 
ral effective models such as the linear sigma model [ll| . 
the Nambu-Jona-Lasinio (NJL) model |12h14| . the chi- 
ral random matrix model [15], chiral perturbation theory 
|16{ and so on. However, those models based on chiral 
symmetry lack any dynamics coming from confinement 
dictated by the Polyakov loop, although there are some 
efforts to clarify the in terplay between chiral dynamics 
and the Polyakov loop |l7l.ll8l|. 

Recent chiral effective models with the Polyakov loop 
degrees of freedom augmented called the Polyakov loop- 
extended NJL (PNJL) model or quark-meson (PQM) 
model J 19. - . 2 6.] are successful from a phenomenological 
point of view to incorporate a coupling between the chi- 
ral condensate and the Polyakov loop. However, these 
PNJL/PQM models are still far from treating the chiral 
condensate and the Polyakov loop on an equal footing, 
except for the work [2^ where the backcoupling of the 
matter sector to the glue sector was discussed by chang- 
ing the phase transition parameter. In fact, the gluonic 
part in these models has several fitting parameters which 
are determined only from lattice QCD data. 

Here we must mention a preceding work for a first- 
principle derivation of confinement/deconfinement and 
chiral-symmetry-breaking/restoration crossover phase 
transition based on the flow equation [23 of the func- 



tional renormalization group |28. I29i given by Braun, 



Haas, Marhauser and Pawlowski Sj] for the full dynam- 
ical QCD with 2 massless flavors (at zero and imaginary 
chemical potential). In this work, the Yang-Mills theory 
is fully coupled to the matter sector by taking into ac- 
count the Polyakov-loop effective potential [30|] obtained 
in a nonperturbative way put forward by Braun, Gies 
and Pawlowski 
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|32| and Marhauser and Pawlowski 
The main purpose of this paper is to provide a the- 
oretical framework (a reformulation of QCD) which en- 
ables one to describe in a unified way the chiral dynam- 
ics and confinement signaled by the Polyakov loop. We 
give an important step towards a first-principle deriva- 
tion of confinement/deconfinement and chiral-symmetry 
breaking/restoration crossover transition. In fact, we 
demonstrate that a low-energy effective theory of QCD 
obtained in simple but non-trivial approximations within 
this framework enables one to treat both transitions si- 
multaneously on equal footing. 

The basic ingredients in this paper are a reformula- 
tion of QCD based on new variables i34 41, 43 46, 48- 
[53j and the flow equation of the Wetterich type in the 
Wilsonian renormalization group |27H29| . The reformu- 
lation was used to confirm quark confinement in pure 
Yang-Mills theory at zero temperature and zero density 
based on a dual superconductor picture [561]. In this pa- 



per, it is extended to QCD at finite temperature and 
density. In principle, our framework can be applied to 
any color gauge group and arbitrary number of flavors. 
For technical reasons, however, we study two color QCD 
with two flavors in this paper. The three color and/or 
three flavor case will be studied in a subsequent paper. 
In future publications, this framework will be applied to 
investigate the QCD phase diagram at flnite density. We 
hope that this paper will give an insight into this issue 
complementary to other works, e.g., J33l [. 

In sec. II, we give a reformulation of QCD written in 
terms of new variables and explain why the reformulated 
QCD is efficient to study the interplay between confine- 
ment and chiral-symmetry breaking. 

In sec. HI, we choose a specific gauge (modified 
Polyakov gauge) to simplify the representation of the 
Polyakov loop. We can choose any gauge to calculate the 
Polyakov loop average and the chiral condensate, since 
both are gauge-invariant quantities and should not de- 
pend on the gauge chosen. 

In sec. IV, we give a definition of the Polyakov loop 
operator and examine how the Polyakov loop average is 
related to the average of the time-component of the gauge 
field. 

In sec. V and sec. VI, we study the confine- 
ment/deconfinement phase transition in pure SU(2) 
Yang-Mills theory at finite temperature. We exploit the 
Wilsonian renormalization group in our framework to ob- 
tain the effective potential VcS of the Polyakov loop L, 
whose minimum gives the Polyakov loo p a verage (L) . It 
is known that the Weiss potential Vw [SOl calculated in 
the perturbation theory to one loop exhibits spontaneous 
center-symmetry breaking, i.e., deconfinement, irrespec- 
tive of the temperature T. This result can be used at 
high temperature where the perturbation theory will be 
trustworthy due to asymptotic freedom, while nonpertur- 
bative approach is necessary to treat the low-temperature 
case. The Weiss potential can be improved according to 
the Wilsonian renormalization group to obtain a nonper- 
turbative effective potential which is valid even at low 
temperature. 

In sec. V, we write down the flow equation of the 
Wetterich type for the effective potential of the Polyakov 
loop in our framework. In fact, the effective potential 
obtained by solving the flow equation in a numerical way 
shows the existence of confinement phase below a certain 
temperature T^. This solution was shown for the first 
time by Marhauser and Pawlowski [3l| and by Braun, 
Gies and Pawlowski [S^], see (SJ] for the previous works. 
In this sense, this section is nothing but the translation 
of their results [3l|, Is^] into our framework. 

In sec. VI, we give a qualitative understanding for 
the confinement/deconfinement transition given in sec. 
V based on the Landau-Ginzburg argument. We answer 
a question why the center-symmetry restoration occurs 
as the temperature is decreased, by observing the flow 
equation for the coefficient of the effective potential. 

In sec. VII, we describe the low-energy effective in- 



teraction among quarks by a nonlocal version of the 
(gauged) NJL model in which the effect of confinement is 
explicitly incorporated through the Polyakov loop depen- 
dent nonlocal interaction. The resulting effective theory 
is regarded as a modified and improved version of nonlo- 
cal Polyakov-loop extended Nambu- Jona-Lasinio (PNJL) 
models proposed recently by Hell, Rossner, Cristofore- 
tti and Weise [22], Sasaki, Friman and Redlich [2J], and 
Blaschke, Buballa, Radzhabov and Volkov [2J|, extend- 
ing the original (local) PNJL model by Fukushima '19'|. 
The nonlocal (gauged) NJL model can be converted to 
the nonlocal (gauged) Yukawa model to be bosonized to 
study the chiral dynamics. 

In sec. VIII, we show that the nonlocal NJL interaction 
among quarks becomes temperature dependent through 
the coupling to the Polyakov loop. This is a first non- 
trivial indication for the entanglement between the chiral 
symmetry breaking and confinement. This feature was 
overlooked in conventional PNJL models. 

In sec. IX, we consider how to understand the en- 
tanglement between confinement and chiral symmetry 
breaking in our framework. This is just a short sketch 
for our strategy following the line given in the preceding 
sections. 

The final section is used to summarize the results and 
give some perspective in the future works. Some technical 
materials are collected in Appendices. 



II. REFORMULATION OF QCD 

To fix the notation, we write the action of QCD in 
terms of the gluon field ^ and the quark field ip: 

•S'qCD =»S'q -I- ^YM , 



^YM :=y d^x— tr(^^,K]^^''[^]), (1) 

where ip is the quark field, ^j = £//^Ta is the gluon 
field with su{Nc) generators Ta for the gauge group G — 
SU{Nc) (A = 1, • • • , dimSU{Nc) = N^ - 1), rhq is the 
quark mass matrix, fiq is the quark chemical potential, 
7'^ are the Dirac gamma matrices (/i = 0, • • • , D — 1), 
Vfj^l^/] := d^ — ig^/^ is the covariant derivative in the 



fundamental representation, ^ 



fiv[ 



:= d„ 



a. 



ig[s!/fj,, s^v\ is the field strength and g is the QCD coupling 
constant. In what follows, we suppress the spinor, color 
and fiavor indices. 

The main purpose of this paper is to give a the- 
oretical framework for extracting a low-energy effec- 
tive theory which enables one to discuss the con- 
finement/deconfinement and chiral-symmetry break- 
ing/restoration (crossover) transition simultaneously on 
an equal footing. We reformulate QCD in terms of new 
variables which are efficient for this purpose. We start 
with decomposing the original SU{N) Yang-Mills field 



£/f,{x) = s:/f{x)TA into two pieces f^, = y^{x)TA and 

<(a;) = r^ix) + ^^{x), (2) 

to rewrite the original QCD action into a new form: 



S. 



d"- 



-i(*g[jr^,jr.])2}+5FP, (3) 

where ^^-^ := g'ipjf^TAiJ is the color current, -D/i['^] := 
dfi — igl^tn ■] is the covariant derivative in the adjoint 
representation and 






/j.f I J ' 



= - {dpS^'' + g/^^^r/)(9''<5^^ + gfCFBypF^ 
- - a^^^^ + g2fAECfBFCyEypF 

+ 2gf^^'^rfdP + gf-^^'^dPyf. (4) 

In what follows we use the notation ^ • ^ for two Lie- 
algebra valued functions ^ = £/^Ta and ^ = ^^Ta 
in the sense that .s^ ■ ^ ■= s/^^^ = 2tr(^^) and 
especially ^^ := ^ • ^Z = si^s^^. 

Historically, the decomposition of Yang-Mills theory 
into new variables has been proposed by Cho [3J| and 
Duan and Ge J35] independently, and readdressed later 
by Faddeev and Nicmi 36] . The decomposition was fur- 
ther developed by Shabanov 37]. 

The decomposition ^ is performed such that 'f^^ 
transforms under the gauge transformation just like the 
original gauge field ^: 

^mW ^ "^''M - ^(^)C^m(2^) + ig-^dp)^-\x), (5) 
while 2^p transforms like an adjoint matter field: 

3C^,{x) -^ SC'^^ix) = ^{x)5i;p{x)^-^{x). (6) 

In the decomposition ([2|), we introduce a new field 

n(x) - n^(a;)T^, (7) 

with a unit length in the sense that n^{x^n^{x^ = 1, 
which we call the color field. In the decomposition ([2]), 
the color field n[x) plays a crucial role as follows. The 
color field is defined by the following property. It must be 
a functional or composite operator of the original Yang- 
Mills field s^p, (x) such that it transforms according to the 
adjoint representation under the gauge transformation: 



n{x) — > n'{x) = fl{x)n{x)n, ^(x). 



(8) 



The color field plays the key role in the reformulation. 
Once a color field is given, the decomposition is uniquely 



determined by solving a set of defining equations and 
hence y^{x) and ^^{x) are written in terms of ^^(x) 
and n{x). For G — SU{2), the defining equations are 
given by 

(I) covariant constantness of color field n{x) in ^^{x): 



Q = D^[r]n{x), 
(II) orthogonality of .^(a;) to n{x): 
0= 3t:^{x) ■ n{x). 



(9) 



(10) 



Then the decomposition for G = SU{2) is uniquely de- 
termined as 



't^ix) = c^{x)n{x) +ig '^[n{x),d^n{x)], 

Cf^{x) := £/f_,{x) ■ n{x), 
^^{x) = tg-^[D^[£/]n{x),n{x)], 



(11) 



To arrive at the result ([3]), we have used the following 
facts. See Appendix [Al for the details. 

(i) The O(^) terms vanish, ^^^"'[y] ■ {D^[r]3C'^ - 
D,y[Y]^f^) = 0, from the property of the new variables 
as shown using the defining equations of the decomposi- 
tion ^. This is somewhat similar to the usual back- 
ground field method in which 0{^) terms in the quan- 
tum fluctuation field .^^ are eliminated by requiring that 
the background field 1^ satisfies the classical Yang-Mills 
equation of motion, i.e., Df^[i^]^'^''[y] = 0. In our case, 
however, 1^ do not necessarily satisfy the classical equa- 
tion of motion. 

(n) To obtain Q^^[y] in dH), an ©(S'^) term is elim- 
inated, -^^^'^D^'^['r]D^^['r]^''^ = 0, by imposing 
the condition: 



£i^[r]jr^ = o. 



(12) 



For the reformulated QCD to be equivalent to the orig- 
inal QCD, we must impose such a constraint to avoid 
mismatch in the independent de gree s of freedom, which 
is called the reduction condition |38l . |41| . 

(iii) The 0(^3) term vanishes, ^{D^[r]^^ - 
D,[y]^^)■^g[3r^',3r'']=0, since D^[r]^,-D4r]^f, 
is orthogonal to [jT'^, Jf'^]. 

For G — SU{2), Y can be chosen in such a way that 
the field strength ■!^\f\ of the field ^ is proportional to 
n: 

^^.m(x) ■.^d,,%{x) - 9.r^(x) - ig{r^[x), %{x)] 

=n{x)Gf^^{x), (13) 

where G^^ is a gauge-invariant antisymmetric ten- 
sor of rank 2, i.e., ^'^,[r]{x) = ^^A^']{x) = 
^ix).^^,^[r]{x)n-^{x) = n'{x)Gf,^{x). The explicit 
form of G|^^^ is written in term of ^^(x) and n(x) as 

Gf,u{x) = df,[n{x) ■ s^i,{x)\ - di,\n{x) ■ ^^(x)] 



In the present approach, we wish to regard the field 
decomposition as a change of variable from the original 
gluon field to new variables describing a reformulated 
Yang-Mills theory in the quantum level [33, 140, Sif (see 
[42-49] for the corresponding lattice gauge formulation). 
To achieve this goal, first of all, n{x) must be written 
as a functional of sip_{x) and thereby all new fields are 
written in terms of the original gluon field J2/^(a;). Such 
a required relationship between s^p{x) and n{x) is given 
by the reduction condition which is given as a variational 
problem of obtaining an absolute minimum of a given 
functional. The condition for local minima is given in 
the form of a differential equation. For G — SU{2), 



[n{x),D^[£/]D^[£/]n{x)]^0. 



(15) 



This is another form of ([T2|) . See 38] in SU{2) case and 
[m in SU{N) case for the full details. 

Remarkable properties of new variables are as follows. 
First, we remind you of the role played by the field Y. 

• The variable 1^ alone is responsible for the Wilson 
loop operator Wc[s^] and the Polyakov loop operator 
L[j2/] in the sense that 

Wc[£/] = Wc[y], L[.s/] = L[r]. (16) 

where the Wilson loop operator is defined by 



Wc[A] :=7V"Hr 



^exp \ig f dx^S!^^{x) 



c 



(17) 



ig ^n{x) ■ [dfj,n{x),d^n{x)]. 



(14) 



where ^ denotes the path ordering and the normaliza- 
tion factor M is the dimension of the representation i?, 
in which the Wilson loop is considered, i.e., J\^ :— dp — 
dim(lij) = tr(l/{). The Polyakov loop operator will be 
defined later. In other words, ^^ do not contribute to 
the Wilson loop and the Polyakov loop in the operator 
level. This is because the defining equation for the de- 
composition is a (necessary and) sufficient condition for 
a gauge-invariant Abelian dominance (or 'f dominance) 
in the operator level. This proposition was first proved 
in 50] for SU{2) and for SU{N) in the continuum [5l| 
and for SU{N) on a lattice (S^]. On the lattice, the 
equality does not exactly hold due to non-zero lattice 
spacing e, but the deviation vanishes in the continuum 
limit of the lattice spacing e going to zero, e — > 0. It 
should be remarked that both the Wilson loop opera- 
tor and the Polyakov loop operator are gauge-invariant 
quantities and that their average do not depend on the 
gauge fixing condition adopted in the calculation. 

• We can introduce a gauge-invariant magnetic 
monopole current k in Yang-Mills theory (without mat- 
ter fields) where k is the {D — 3)-form. For Z? = 4 and 
G = S'C/(2), 

k^{x) :-a,*G^,(a:), (18) 

Gpi. :=n • ^^,y[y] = d^c^ - 5^c^, + ig^^n ■ [S^n, di,n], 

where /^^ is gauge-invariant field strength. This is 
because the field strength ^^,^[i^] := d^% - d^^Y^ - 



*ff[^i ^] is proportional to n: 

■^ii^v\y\ = n{df,c^ - d„Cf, + ig^^n ■ [df,n,d^n]}. (19) 

• The gauge-invariant "Abelian" dominance (or 'f dom- 
inance) and magnetic monopole dominance (constructed 
from "f) in quark confinement liave been confirmed at 
T = (and /ig = 0) by comparing string tensions cal- 
culated from the Wilson loop average by numerical sim- 
ulations by [3 for SU(2) and by [^ for SU(3). Here 
it should be remarked that the "Abelian" dominance is 
the dominance for the vacuum expectation value (or av- 
erage) : 



(PFcKl) ^ {Wc[y]), (iH) c, {L[r]). 



(20) 



Next, we pay attention to the role played by the re- 
maining field ^ . 

• In the absence of dynamical quarks (corresponding 
to the limit niq = oo of QCD, i.e., gluo dynamics), 
^if decouple in the low-energy regime as the correlator 

{^if'{^)'^if{y)) behaves like a massive correlator with 
mass Mx- In fact, numerical simulations demonstrate 
for G = SU{2) and L> = 4 m 



Mx = 1.2 - 1.3 GeV. 



(21) 



We can understand this result as follows. The field ^^ 
can acquire the (gauge-invariant) mass dynamically. This 
comes from a fact that, in sharp contrast to the field ^^, 
a "gauge-invariant mass term" for .^^ can be introduced 



1 



M 



^^ (x)^^ [X), 



(22) 



since ^.f 



{x).^if-{x) is a gauge-invariant operator. More- 
over, this mass term can originate from a vacuum con- 
densation of "mass dimension-2" , (^^^ {x) ^^ (x)') ^ 
as proposed in [55|. In fact, this condensation can be 
generated through self- interactions 0(^'*) among J^ ,. 
gluons, M| ~ {^J^{x)^f{x)), as examined in ©lij 
It is instructive to remark that the value ([21]) agrees with 
the earlier result of the off-diagonal "gluon mass" Ma in 
the Maximally Abehan (MA) gauge ^SCi] for 51/(2) case. 
Ma ^ 1.2 GeV. See ^ for SU{3) case. Ma ^ 1.1 GeV. 
In MA gauge, it was shown that even at finite tempera- 
tures Abelian dominance (diagonal part dominance and 
off-diagonal part suppression) holds for the spatial prop- 
agation of gluons in the long distance greater than 0.4fm. 
It was observed that the diagonal gluon correlator largely 
changes between the confinement and the deconfinement 
phase, while the off-diagonal gluon correlator is almost 
the same even in the deconfinement phase [62| . Although 
the similar results are expected to hold in our formula- 
tion, this observation must be checked directly, as will be 
confirmed in |63| . 

• In the presence of dynamical quarks {niq < oo), ,^ 
is responsible for chiral-symmetry breaking in the follow- 
ing sense. We consider to integrate out the field ^"^ in 
a naive way. This helps us to obtain an intuitive and 



qualitative understanding for the interplay between the 
chiral symmetry breaking and confinement. Later, this 
integration procedure will be reconsidered from the view- 
point of the renormalization group to obtain a systematic 
improvement of the result. 

Here we neglect 0{^^) and 0{^'^) terms, which will 
be taken into account later. Then the integration over 

^,f can be achieved by the Gaussian integration ac- 

^ ^ — ii 

cording to :,4Q]. Consequently, a nonlocal 4 fermion- 

interaction is generated: 



^QCD 



-S"' 
-'-', 



glue 
off ■ 



Si 



gNJL 
off ' 



J, 



glue 
eff 



= / d^x—{^,An)' 



In det Q [r] ;^^^ - i In dot G [r ] -^^ , 



S: 



■gNJL 
off 



d^x {pii-ft'V^ir] -mg + i7V)V' 



d^a 



(23) 



where the last term — In det G in S^g° comes from 
the Faddeev-Popov determinant associated with the re- 
duction condition ([T^ (see ^j for the precise form). 

This is a nonlocal version of a gauged NJL model (re- 
alized after Fierz transformation). The chiral-symmetry 
breaking/restoration transition and the phase structure 
of a local version of gauge NJL models were first stud- 
ied by solving the Schwinger-Dyson equation in the lad- 
der approximation for QED-like |64| - [6a | (see 67 1 for a 
review) and QCD-like [68] running gauge coupling con- 
stant. They are confirmed later by a systematic approach 
of the renormalization group [69| . 

The range of the nonlocality is determined by the cor- 
relation length ^, which is characteristic of the color ex- 
change through gluon fields. Therefore, this correlation 
length ^ is identified with the inverse of the effective mass 
Mx,i.e.,^~M ' 
correlator, see (jSl) 

In other words, Mx is identified with the ultraviolet 
cutoff A below which the effective NJL model appears 
and works well. Interesting enough, Mx is nearly equal 
to the ultraviolet cutoff adopted in the NJL model 

y^< A4 = 1.4GeV, IpI < A3 = 0.6GeV, (24) 

see [13. 

We can decompose the gauge field jz/^ into the the low- 
energy (light) mode "p < Mx and high-energy (heavy) 
mode p > Mx ■ 

- p") ^ j^,{x)e{p' 



^. In fact, {Q-^)'^^{x, y) is the ^ field 



Kip) =<(p)e(M 



-Mi; 



(25) 



In the above treatment, ^fj,(j)) is supposed to have only 
the high-energy mode. The low-energy mode, if any, will 
be responsible for the vacuum condensation [40'|. For 
the precise understanding, we need the renormalization 
group treatment as given later and the implications for 
the nonlocal NJL model will be discussed there. 



III. GLUON SECTOR AND GAUGE FIXING 

The Polyakov loop operator L and the chiral operator 
ip^l) are gauge-invariant quantities. Therefore, their aver- 
age do not depend on the gauge-fixing procedure adopted 
in the calculation. We can choose a gauge in which the 
actual calculation becomes easier than other gauges. 

In what follows, we treat the time-component Yq 
and space-component Yj of Y^ differently to consider 
the finite-temperature case. We consider the following 
Polyakov gauge modified for new variables in our refor- 
mulation. If the color field n-^{x) is uniform in time, 

aon^ (s) = ^ n-^ {x) ^n'^ix), (26) 

then Yq reduces to 

Y^''{x)^c^{x)n^{x) (A = 1,2, 3). (27) 

Moreover, if co(x) is uniform in time, 

(9oCo(a;) =0^co(a;) =co(a?), (28) 

then Yq reduces to 

Y^''{x)^co[x)n^{x) (A = 1,2, 3), (29) 

which satisfies 

aoV(a;)=0 (A = 1,2, 3). (30) 

In this setting, Yj are given by 

Yj'^ix) = Cj{x)n'^{x)+g-h^^^djn^{x)n^{x). (31) 

(1) In order to simplify the calculation of the Polyakov 
line, we adopt the Polyakov gauge in which the gauge 
field is diagonal and time-independent: for the back- 
ground field ^^(a;), 

V(a;)=co(a;)<5^^ 

which leads to 

ao^o^(x) = o. 

This is realized, if we take the gauge ^ 

n^{x) = S^^ (34) 

In this gauge, the space-component reads 



(32) 
(33) 



Y/'ix)^c,ix)S-'. 



(35) 



^ This is an oversimplified choice for the color field 'n{x). By 
this choice, we can not separate the non-perturbative contribu- 
tion coming from topological configurations such as magnetic 
monopole. It is desirable to take into account color field degrees 
of freedom explicitly to see the effect of magnetic monopole in 
the confinement/deconfinement transition. 



which is not time-independent, doYj^{x) ^ 0. 

(2) We expand the theory around the non-trivial uni- 
form background g^^TipS"^^ for the time-component Yn^, 
while the trivial background for space-components Y, -."^ 

V(x) ^co{x)S^^ co{x)=g-'Tip + vo{x), 
Y^{x)=0 + vf{x), (36) 

such that {co{x)) = g^^T{ip) + {vn{x)) = g^^T{ip) with 
{vo{x)) = and {vl^{x)) = 0. Here the prefactor g^^T = 
{gf3)~^ was introduced just for the purpose of simplifying 
the expression of the Polyakov loop, see (|^T|) . 

(3) We take into account the expansion up to quadratic 
in the fluctuation fields Vq and Vj, which we call the 
quadratic approximation. 

In the calculation of Q^^, if we neglect all fluctua- 
tion fields vq and Vj, namely, Y^{x) — g^^T{pS^^ and 
Y/-{x) = 0, then we can put ^^^['^l = in Q^^ , and 
Qfjf is diagonal in the Lorentz indices: 

Qt."" = G^'^a.u, (37) 

In this approximation, we have 

qAB ^ ^sABq2 _ jjACp^^j^CB^^^^ (38) 

with 

~D^'^[Y]D^''[Y] = ~ d^'^dl + 2e^''^T^do 

+ {5^B -5^H''^){T^f. (39) 

Thus we rewrite the gluon part 5'^g"''[l^] as 



tr\ 



1 



(i / d?xYo{x){-djdj)Yo{x) 



+ i j d^xY^{x){^8^''dj - S^^'dDYrPix) 
+ \ I d^^Y^^{x){-5^^dl}Y^{x), 



AB 



- In det Q[r];^j' - i In det Gir]-^"", 



(40) 



where /3 is the inverse temperature /3 := l/T, and Yt and 
Yl denote the transverse and longitudinal components of 
T^i respectively. 



^ I have assumed that the spatial component of ^ has a triv- 
ial background. In view of logical consistency, one must expand 
the spatial and temporal components around non-trivial back- 
grounds, and then one must search for the minima of the effective 
potential calculated as a function of two variables, i.e., the tem- 
poral and spatial backgrounds. In this paper, it is assumed that 
a minimum is realized at vanishing spatial background and that 
the neglection of the spatial background does not so much affect 
the confinement/deconfinement transition temperature. Indeed, 
it must be checked whether this assumption is good or not. 



IV. POLYAKOV LOOP 



For G = SU{2), the Polyakov loop operator L{x) 
L['fQ{x, •)] is defined by 



Lix) :=-tr(P), 
P{x) -.—^exp ig 



/3 = 1/T 



TA 



da;oro^(a;,xo)^ 



(41) 



where PP^ = 1 and detP =1. In the above gauge 
choice, 



P(a;)=exp ig(3co{x)n (x) — 



After a suitable (i-independent) gauge transformation, 
the color field n'^{x) is eliminated: 



L{x) = -tr(exp igPcoix)-^ ) = cos f 1 



(42) 

Owing to periodicity and center symmetry, we can 
restrict the Polyakov loop average to (L) > for 
G = SU{2). Then the Polyakov loop average {L[r]) 
is bounded from above by L[('^(a;, •))]: 



0<{L[%{x,-)])<L[{%ix,-))] 



{^) 



(43) 



where we have only to consider the range < ip < tt. 
This inequality follows from the Jensen inequality, since 
cos(a;) is concave for < x < 7r/2, see |3l| . 

In the case of ruq = oo , if the center-symmetry is bro- 
ken (L) > 0, namely, deconfinement takes place, then the 
vacuum (as a minimum of the effective potential T4ff ((p)) 
is realized at {ip) < tt. If the vacuum is realized at 
(if) = TT, then the center-symmetry is restored (L) = 0, 
namely, confinement occurs. The relation (j42p yields the 
relationship for the average between the gauge field and 
the Polyakov loop operator: 



(arccosL(a;)) = 



gPMx)) 



2 



^, (44) 



where the left-hand side is the average of an gauge- 
invariant object (since L is gauge invariant) and happens 
to agree with the average {'Vq) of the gauge field in the 
Polyakov gauge. It is also shown [31[ that the converse 
is true: In the center-symmetry-restored phase, (ip) = tt, 
since 

= (arccosi(a;)) = arccos(i(a;)) = — . (45) 

Therefore, (Yq) or (ip) in the Polyakov gauge gives 
a direct physical interpretation as an order parame- 
ter for the confinement/deconfinement (order-disorder) 
phase transition. The effective potential Ucs{{L)) of the 
Polyakov loop average (L) could be different from the 



effective potential Vcff(('^)) of the gauge field average 
{'fo) in the following sense. Although both potentials 
give the same critical temperature Td as a boundary be- 
tween (L) = and (L) ^ 0, the value of the effective 
potential Ucs{{L[%])) does not necessarily agree with 
UcsiLKYo)]) = Kfr((%)) at a given temperature T, since 
we have only an inequality {Llfo]) < L[(7^)]. This differ- 
ence could affect the critical exponent and other physical 
quantities of interest. Therefore, the result obtained from 
ycs{{%)) must be carefully examined. 



V. DERIVING THE 

CONFINEMENT/DECONFINEMENT 
TRANSITION 

In this section, we restrict our consideration to the pure 
glue case. We show that the pure gluon part S'fA° can de- 
scribe confinement/deconfinement transition signaled by 
the Polyakov loop average (L). In this section, wc com- 
pletely follow two remarkable papers by Marhauser and 
Pawlowski 31] and by Braun, Gies and Pawlowski [32|, 
which succeeded to show the transition for the first time 
based on the functional renormalization group (FRG). In 
the next section, we explain how these results are under- 
stood from the Landau- Ginzburg argument. 

We consider the flow equation called the Wetterich 
equation [23] for the fc(RG scale)-dependent effective ac- 
tion r^: 



dtT^m =-STr 






nm 






i?$ , 



dtR'S>.. 



(46) 



where t is the RG time t := In -^j 9t := ^ — k-^ ior some 
reference scale (UV cutoff) A and R<s>^k is the regulator 
function for the field $. Here STr denotes the super-trace 
introduced to include both commutative field (gluon) and 
anticommutative field (quark, ghost). See [2^ [23| for 
reviews of the functional renormalization group. 

If we restrict our consideration to the pure glue case 
^YM under the gauge n"^{x) = 6^^, then the relevant 
fields $ are y^^ix) , ^i^{x) and FP ghosts (ghost and 



antighost) '^^{x),'i^{x), i.e., $t = {yA^ ^a^c^a^c^a 
In this section, we use the Euclidean formulation 



In 

the modified Polyakov gauge and within the quadratic 
approximation adopted in sec. IIIIl 



dtV 



Ufc 



-Tr 



ri fe- 



H -1 



(5rt "^sr 



Rk 



■dtRu 



dt-Tr{HQ^^ 



S^^S,, 



9tTr{lnJG 



AB 



■fj,iyRk]} 

S^'^Rk]}. 



(47) 



where the second contribution in the right-hand side 
comes from the ^ field and the last one from the ghosts 



fields 39| , and we have used the same regulator function 
Rk for the gluon and ghost up to the difference due to 
the tensor structure. 

We neglect back-reactions of the Yq potential on the 
other gauge fields Yj , as in the treatment |31i] . Assuming 

an expansion around 'fj = 0, F^ ^ :— j^^k-^ is block- 
diagonal like the regulators, and the flow equation can 
be decomposed into a sum of two contributions: under 
the approximation ([57]) . 



(normalized) SU(2) Weiss potential V„.(^) 



dtTk =^Tt 



1 



.(2) 



■dtR 



rr + Rk 

atTr{ln[G"^^ + S^^ 



't-nk,fj,iy 



Rk]}, 



(48) 



where the gluon regulator Rk.^iv is a block-diagonal ma- 
trix in field space, 

-Rfe,00 =Ro,k = ZoRopt,k{P ), Rkfij = = RkJQ, 

Rk,ji ^RT,kTji[p) = Z,T,,(p)i?opt^kT(p'), (49) 

where Tji := 5jg — ^^ is the transverse projection op- 
erator and i?opt k(p^) is the (3 dim.) optimized choice 



i?opt,k(p')-(fc'-p')0(fc2-p'). 



(50) 



The first term in the right-hand side encodes the quan- 
tum fluctuations of %, while the second one encodes 
those of the other components of the gauge field and 
ghosts. In the present truncation, the second term is 
a total derivative with respect to t, and does not re- 
ceive contributions from the first term. Therefore, we 
can evaluate the flow of the second contribution, and use 
its output Vr^ki'^o) as an input for the remaining flow. 



^.r 



1 



'ti-k 



-J 



(Pp 



1 



r 



(2) 



Rr 



dtRo,. 



(51) 



2" 7 (27r)3 
+ dtVT,k, 
where for uj — 2T:Tn 

-S^^'Rk]} 
rtrln[G^^(w,p) + S^^ {k^ - P^)0{k^ 



Vr.k ■■= TrjlniG-^^ 
(fp 



=^E 




FIG. 1: (The upper panel) (normalized) SU(2) Weiss poten- 
tial Vw as a function of ip. (The lower panel) SU(2) Polyakov 
loop 1/ as a function oi (p, L = cos (^). 



was obtained by one- loop calculation [30J: 



VM/=Trln[G 



AB^ 



=^E 

nez 



d^p 

_fp_ 

(Pp 



(2^) 



trln[G^^(po = c^,p)] 

trln[p2 + (a; + T(^)2] 



(53) 



where we have neglected the (^-independent (or "f^ inde- 
pendent) contributions. 

The closed form of the Weiss potential is obtained after 
summing up the Matsubara frequencies: 



(27r)3 
dPp 

.k^ 



2\i (mod 27r). 



-\{^--f 



1 ( .A . I" 

127r2^^ ' 12 



(54) 



trln[,5^^p2 _ ^2 ^ ^AB(^2^ _ p')0(4 

^E4-r^trln[^-'^^ 
-T^47r 



rC'Y' -^n 



nGZ 



dpp^ 



trln[J'^V-£'o] + Hv. (52) 



Here we have introduced the Weiss potential Vw which 



The Weiss potential Vw is g^ independent and the over- 
pm curve scales as T"*. Vw has symmetries: Vwi—f) — 
Vw if) and Vw {f + 27rn) = Vw (<p) ■ Vw {v) lias minima 
at (/3 = 27rn, and the Polyakov loop has the nonvanish- 
ing value L = cos ^ = (—1)", implying deconfinement. 
See Fig.[T] Therefore, Vw{v) is considered to be vahd at 
very high temperature where the perturbation theory is 
trustworthy. In Fig. [21 we observe 

hmVr.fc = Vt,o = Vwi lim Vr.k = 0. (55) 

k^O ' fc-f-oo 
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FIG. 2: Vr^fc for different values of fc [reprinted from [3i| 



p [GeV] 



After integrating over the fields other than "%, we are 
lead to the effective action of 1^, 



^k\M =/? / d'^ \ -\za%{x)d,d,%{x) + v^^ 



v!tk[n^vTA%] + ^yk[%] 



(56) 



Then the flow equation is reformulated for AVfe with the 
external input VtX- 



mAv,m = lpJ^ 



1 



.(2) 



Rk, 



dtRi 



't^O.k 



00 



where 



r['>[%]^f3{z,p' + d\Vk[%]} 



(57) 



(58) 



Using the specific regulator, i?o,fe = Zo{k'^ —p'^)6{k'^ — 
p^), which yields 

dtRo.k = [dtZoik^ - p') + 2Zofc2] e{k'' ~ p"), (59) 

we can perform the momentum integration analytically. 



pdt(^Vk[%]) 
2 1 



(?7fc/5 + l)fc5 



"3 (27r)2 Z-^g^fi^dl{VT,km + ^Vk[%]) + k^ 



, (60) 



where we have introduced the running coupling ak de- 
fined by 



9l 



.-i9 



9i--=^k'9'. «fc:=^ = ^/f^, (61) 

and the anomalous dimension -qk defined by 

77fc ::= dt InZfe = -dt Ina^. (62) 

By introducing the dimensionless RG scale k and the 
dimensionless effective potential V defined by 



FIG. 3: The running gauge coupling constant Os for temper- 
atures T = 0, 150, 300, 600 MeV [reprinted from [sT 



-0.1 - 



<> 



-0.2 - 




T = 500 MeV 

T = 400 MeV 
T = 350 MeV - 
T = 300 MeV 
T = 250 MeV 



■n.12 



3/2 7U 



271 



FIG. 4: Full effective potential VJ,|f"°, normalized to at (/p = 
[reprinted from [3ll|]. 



the flow equation is simplifled as 



dt^Vd%] 



(1 + r;fc/5)fc2 



67r^ 1 



fc2 "¥> 



di{v^t[%]+^Vi^[%]y 



(64) 

where all scales are measured in units of temperature. 
It turns out that the input in solving the flow equation 
is just a running gauge coupling constant au- A spe- 
ciflc choice for the running gauge coupling constant is 
given in Fig. [3] For the derivation from the renormaliza- 
tion group, see |71| . The flow is initialized in the broken 
phase at any temperature. By solving the flow equation 
in a numerical way with an input for the running gauge 
coupling given in Fig. |3l the full effective potential T4ff 
(normalized to at (/? = 0) is obtained in Fig. |3] for vari- 
ous temperature. 

According to f3l|, a second order phase transition oc- 
curs at a critical temperature 



k:=(ik^ k/T, V := I3^V = V/T^ 



(63) 



T, 



305;^°MeV, Td/V^ 



0.69; 



-.04 
-.12' 



(65) 
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with the string tension ^/a = 440 MeV. This agrees 
within errors with the lattice resuh Tdl\fo — 0.709. 
Moreover, these results were confirmed by considering 
another gauge |33 |. 



VI. UNDERSTANDING THE EXISTENCE OF 

CONFINEMENT TRANSITION ACCORDING TO 

THE LANDAU-GINZBURG ARGUMENT 

In this section, we show that some qualitative aspects 
of the deconfinement/confinement transition found in the 
previous section can be understood without detailed nu- 
merical works, although the precise value of the transi- 
tion temperature T^ cannot be determined without them. 

For G — SU{2) in the pure Yang-Mills limit mq — > 
cx), the effective potential Vg\-ac{L) for the Polyakov 
loop L must be invariant under the center symmetry 
Z(2). Therefore, Vgiuc{L) is an even function of L, i.e., 
Vg\nc[L) — Vgiue(— i) where L is real- valued L = L* . 
Thus the Landau-Ginzburg argument suggests that the 
effective potential V^fj™(L) for G = SU{2) has the power- 
series expansion in L near the transition point L = 0: 



^^t' 



vr^'iL) 



C + ^L^ + ^L^ + OiL'). (66) 



As the vacuum is specified as minima of the effective po- 
tential, the conflnement/deconfinement transition tem- 
perature Td is determined from the condition C2{Td) ~ 
so that the low-temperature [T < Td) confinement phase 
(L) = is realized for C2(T) > 0, while the high- 
temperature (T > Td) deconfinement phase (L) ^ 
is reahzed for C2(T) < 0, provided that the positivity 
Ci{T) > is maintained across the transition tempera- 
ture. Consequently, the transition is of the 2nd order. 

Indeed, we confirm that the Landau-Ginzburg 
description is correct and valid for the confine- 
ment/deconfinement transition, by making use of the 
fiow equation given in the previous section. This is a 
microscopic justification of the Landau-Ginzburg argu- 
ment for the confinement/deconfinement transition. In 
our treatment, however, it is more convenient to write 
the effective potential V^fj"" in terms of the angle vari- 
able (p (rather than L) around the transition point ip — t: 
(instead of i = 0). Defining (^ := (p — tt, we find that 
Vjj^''(i^) must be an even function V^^°{(p) = V^f^'^{—(p) 
due to the center symmetry and hence odd terms (e.g., 
(p, iffl) do not appear: 

y,r(^)=co+^95^. 






0{^% 



(67) 



At sufficiently high temperature 

ECIUC 



we observe that 



C2{T) < and hence V^g"° has the minimum at ip y^ 
(-4=> i 7^ 0) leading to deconfinement. In order to 
show the existence of the confinement/deconfinement 
transition at T = Td, C2{T) must change the signature 
C2{T) > below this temperature T < Td and hence 
the minimum occurs at ip = (<^=^ L = 0) leading to 



confinement. Therefore, the confinement/deconfinement 
temperature Td is determined by G2{Td) — 0, provided 
that the positivity C4(T) > is maintained. 

For this purpose, we study the scale dependent effec- 
tive potential V^g^^ at fc > 



V: 



glue _ 



cff.k 



„ G2,k -2 , ^4,k _4 „ , ,i 



^6\ 



(68) 



2 '' ' 4! 

and see how it evolves towards the limit fc — > accord- 
ing to the fiow equation to obtain the physical effective 
potential F/^^ -Jet^^o- 

As in ([5B)). Vjf™ is decomposed into two pieces: 



V: 



glue 



off,k 



Vn 



T,k 



^v,, 



(69) 



where we have defined the dimensionless potential ac- 
cording to the rescaling ([55)) . The first part V^^ is the 
(fc-dependent) perturbative part (|52|) obtained essentially 
by the one-loop calculation with the regulator function 
Rk being included. For this part, the closed analytical 
form can be obtained, see Appendix IB] While the second 
part AVj^ represents the non-pertubative part which is 

initially zero AVj,|fc=A = and is generated in the evolu- 
tion of the renormalization group. This part is obtained 
only by solving the fiow equation (j64p and its analytical 
form is not available (at this moment). 
We expand Ij, ^ in powers ot (p — p ~ n: 



V. 



T,k 



A I "^2,fc -2 



^4,fe -4 



o(^'), 



(70) 



where coefficients are drawn as functions of k in Fig. [SJ 
see Appendix [Bj for their closed analytical forms. 

A coefficient A^- of ip^ in V r 
^10 12 



A coefficient A ■ of ip in V - 




FIG. 5: A2,k and A4^k as functions of k. 

Suppose that AVj. is of the form: 

AV-,^ao,k + ^p' + ^^' + Om. (71) 
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A flow equation ((64|) for the effective potential ((56|) is 
reduced to a set of coupled flow equations for coefficients 
in the effective potential ^ with ([70|) and (fTTI) : 



anomalous dimension 



9^02, fe 



9^04, fe — 



(1 + 


i^.)fc2 


47rafc 

fc2 


(A4., 


+ a4^fe) 


67r2 


[1 + ^ 
Gr47rQfc 


^(A2,fc + a2,fc)]2 
(A4^fe+a4^fc)]2 



6772 



[l + ^(A2,fc + a2,fc)P 



(72) 



which are coupled first-order ordinary but nonlinear dif- 
ferential equations for coefficients. In Appendix [C] we 
see that this form (|7T|) is justified as a solution of the 
flow equation. In fact, it is easy to see that 9j.ai.fc = 
and dj^as^k — are guaranteed from the flow equation, 
if the effective potential has no odd terms at arbitrary 
k. Therefore, if an initial condition, oi.fe = = 03, ^ 
at fc = A is imposed, then ai^k = and as.fc = are 
maintained for < A; < A by solving the flow equation. 
In performing numerical calculations, however, one must 
truncate the infinite series of differential equations up to 
some finite order to obtain manageable set of equations. 

The running gauge coupling constant 




0.0 



0.5 1.0 1.5 2.0 2.5 

Tlie running gauge coupling constant 




FIG. 6; The running gauge coupling constant a^ at T = 
0.001,0.125,0.25,0.50,1.0 GeV, from the top at the lowest 
temperature T = O.OOlGeV to the bottom at the highest 
temperature T = 1.0 GeV. (The upper panel) Qfe as functions 
of k. (The lower panel) a^ as functions of k. For a given 
temperature, there is a critical value fee separating the deep 
IR region (|76p from the higher momentum region (|75|) . The 
discontinuity of the derivative seen at kc comes from a crude 
approximation in which we have taken into account just the 
first linear term (i.e., ci = C2 = • • • = 0) in the expansion (|76|l . 
and can be avoided if we take into account higher order terms 
as explained below (|76|l . However, this is not essential to see 
qualitative behaviors of the solution of the flow equation. 
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FIG. 7: The anomalous dimension rjk as functions of k at 
T = 0.001, 0.125, 0.25, 0.50, 1.0 GeV. In each graph for a given 
temperature, there is a critical value kc of k separating the 
deep IR region rjk — —1 from the higher momentum (inter- 
mediate and UV) region rjk > 0. The temperature is distin- 
guished by kc ranging from the smallest value at the highest 
temperature T = l.OGeV to the largest value at the lowest 
temperature T = 0.001 GeV where 77^ ~ —1 for k < kc and 
rjk — for k > kc- The discontinuity of the derivative seen at 
kc is due to the same reason as that explained in Fig. [S] and 
is not essential to see qualitative behaviors of the solution of 
the flow equation. 



We can understand qualitatively why a 2nd order 
phase transition from the deconfinement phase to the 
confinement phase can occur by lowering the tempera- 
ture. 

The flow starts from a2.fc = and hence C2,k = 
^2,fc + a2,fe < (because of ^2,*; < 0) at fc = A >• 1. 
We assume C4,fc = A4^k + cn.k > for < fc < A, as a 
necessary condition for realizing a 2nd order transition. 
Otherwise, we must consider the higher-order terms, e.g. 
0{lp^). (This assumption is assured to be true by numer- 
ical calculations of the full effective potential [31, 32], as 
reproduced in the previous section.) This assumption 
allows us to analyze just one differential equation for ob- 
taining qualitative understanding: 



9^02, fc = 



(1 + I^fc) 



6772 [1 



47rafc(A4^fc H- 04, fc) 
+ ^{A2,k + a2,k)? 



(73) 



Then the right-hand side of ([75)) is negative, since the 
running coupling constant ctk is positive and 1 -I- g?7fc 
is positive, see Fig. |6] and Fig. [T] Consequently, 02, fc 
started at zero becomes positive a2.fc > just below A 
and increases (monotonically) as k decreases. See Fig. HI 



Note that the denominator can vanish 1 



-{A 



2,fc ' 



a2,fe) = at some k* (since C2,fc = A2,fc -I- a2,fc < or 
< a2,fc < — ^2,fe) where the right-hand side of (|73)) 
becomes negative infinity and 02, fc blows up there. To 
avoid this pathology and to obtain the solution all the 
way down to the limit fc — > 0, 02, fc must grow relatively 
rapidly so that |^2,fc + a2,fc| ^ 1 towards the limit fc — >■ 0. 
An important observation of the flow equation (l73l) is 
that the explicit temperature-dependence comes from the 
running coupling constant alone. At zero temperature, 
the running coupling constant is well parameterized by 
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the fitting function [72 



ak 



^g A coefficient a^. of (p^ in AV^ ^ 



47r X 0.709/A^c 



In[e + ai(fc2)a2 +6i(fc2)''2 



(74) 



where ai = 5.292, 02 = 2.324, 61 = 0.034, 63 = 3.169. 
in units of GeV. 

For the perturbative region k ^ T, i.e., fc ^ 1, we 
adopt this form: k^ —T^k^, 



«. = # 



An X 0.709/A^c 



47r ln[e + ai(T'^k^Y^ + 6i(T2fc2)b2] 



(75) 



For the nonperturbative region k < 2ttT, i.e., fc < 0(1), 
we adopt the running coupling which is governed by an 
infrared fixed point |71j 



ak —a 



3d 



Cl 



=a 



3d 



k_ 
k + cik^ + C2P 



C2 



(76) 



where coefficients ci , C2 , ... are determined such that the 
coupling at zero temperature (j75p and its derivative with 
respect to k are connected continuously with this ansatz 
([75| at the scale set by the lowest non-vanishing bosonic 
Matsubara-mode uj — 2t:T. However, the running cou- 
pling constant at small momenta ([75| does not contribute 
to the explicit T-dependence in the scaled flow equa- 
tion, since it is written in terms of the scaled k alone 
and hence denoted by a common curve going through 
the origin for any temperature T in the second figure of 
Fig. [Bl Therefore, the running coupling at very small 
momentum region can not be responsible for the confine- 
ment/deconfinement transition at finite temperature, if 
this observation is correct. As can be seen from the sec- 
ond figure of Fig. [HI the dominant contribution comes 
from the intermediate momentum region above 0(1) 
GeV. Thus, we can avoid the issue of gauge-fixing ar- 
tifact in the deep IR region due to Gribov copies in the 
zero-temperature case, see e.g. [Tj, [TJ] and reference 
therein. 

We consider a solution of the reduced (or normalized) 
flow equation as a function of k, rather than fc, for a given 
temperature T. Then the difference between high and 
low temperature phases attributes to the behavior of the 
running coupling constant ak as a function of fc = Tk, 
which brings the explicit T dependence to the reduced 
flow equation. In the case of high-temperature T ^ 1, 
k — Tk becomes large for a wide range of k and the 
running coupling constant ak remains relatively small. 
The resulting slow increase of 02, fe keeps a2,k small such 
that C2^k = ^2,*: + 02, fc < or a2,k < "^2,fe even at 
k = 0. This leads to the center symmetry breaking at 
high-tcmpcrature. 

In the case of low temperature T <C 1, A: = Tk be- 
comes small for the same range of k and the running 
coupling constant ak gets into the intermediate region 
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FIG. 8: 02 k vs. —A2 k as functions of k for T < Td- 



of 0(1) GeV rapidly and becomes larger as the temper- 
ature becomes smaller. At sufficiently low temperature, 
02, fc increases in decreasing k so rapidly that 02, fc eventu- 
ally reaches to the point A2^k + 02, fc = or 02, fc — — ^2,fc 
at a certain value k = ko. In other words, the graph of 
02, fc intersects with that of —A2^k at k — ko. In the region 
< A: < fco where C2,fc = ^2,fc + a2,fe > or a2,fc > -A2.k, 
the flow equation reads 



dta2^k 



(1 + |?7fc) (^4,fc + a4,fc)fc'^ 
67r2 47rafe(A2,fc -|-a2,fc)2' 



(77) 



the right-hand-side gets small negative, and 02, fc becomes 
flat near the IR limit. See Fig.[S] Finally, 02, fc reaches the 
value reahzing C2,fc = ^2,fc + 02, fc > or 02, fc > — ^2,fc at 
k — 0. This leads to the recovery of the center symmetry. 
The difference is clearly seen from the second figure of 
Fig. [S] where the running gauge coupling ak is drawn as 
a function of k for various temperatures. 

In our treatment, the difference between the three- 
dimensional RG scale kx and the four-dimensional one 
k is neglected by equating two scales kx — k just for sim- 
plifying the analysis, since it is enough for obtaining a 
qualitative understanding for the transition. This is not 
be the case for obtaining quantitative results, see Ap- 
pendix C of [31| for the precise treatment on this issue. 



VII. QUARK PART AND GAUGED 
NONLOCAL NJL MODEL 



We examine the quark self-interaction part Sh 
Sd^xSd''y\/^\x)gHQ-^[n)t^{x,y)/'^^{y). 



estimating the effect of Q 
proximation as the above. 



In 

[f], we take the same ap- 
Consequently, the inverse 
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{Q ^)f^J^['f] is diagonal in the Lorentz indices: 

^^m1[F^-F_^] \[F., + F^.,] 1, (78) 
Fo, 



where F^:, is defined by 



F^{id) 



1 



1 



(^a^)2 + (T<^)2 + 2T(^^9o■ 



(79) 



In what follows, we consider only the diagonal parts of 
{G~^)^^ [Y]. This is achieved by the procedure 



9 



{Q-X^^{x,y)^g'^''6^^g{x-y), 



which yields 

Then the nonlocal interaction is obtained as 



G{x - y) = 



2 m^ 



g^ F^ + F_ 



1 



(80) 



(81) 



(82) 



This approximation is used just for simplifying the Fierz 
transformation performed below and hence it can be im- 
proved by taking into account the ofi-diagonal parts of 
G~^ if it is necessary to do so. 

For Z? = 4, we use the Fierz identity [73 to rewrite the 
nonlocal current-current interaction as 



Sint ^ I d^xj d^y/'^^{x)g{x ~ y)^^^{y) 

d'^x / d*yg{x - y)^Caiip{x)raipiy)) 

-' a 

X i,p{y)T^i;{x)) 
= I d'^x I d^zg{z) Y^ Ca{ij{x + z/2)Tc,i^{x - z/2)} 

x{i!{x-z/2)T^i!{x + z/2)}, (83) 

where the Fq, are a set of Dirac spinor, color and fla- 
vor matrices, resulting from the Fierz transform, with 
the property 7oFJj7o — F^. Although the Fierz trans- 
formation induces mixings and recombinations among 
operators, the resulting theory must maintain the sym- 
metries of the original QCD Lagrangian. A minimal 
subset of operators satisfying the global chiral symme- 
try SU{2)l X SU{2)ji which governs low-energy QCD 
with two-flavors is the color-singlet of scalar-isoscalar and 
pseudoscalar-isovector operators Thus, by restricting Fq, 
hereafter to 



and ignoring other less relevant operators (vector and 
axial- vector terms in color singlet and color octet chan- 
nels), we arrive at a nonlocal gauged NJL model 

5f^JL ^ j d^x4'{x){iY'V^,[r] -rn^ + il^^mx) + 5i„t, 

S'int = f d^x f d^zg{z)[{p{x + z/2)F„V(a; - z/2) 

xi;{x- z/2)rc-(p{x + z/2)]. (85) 

This form is regarded as a gauged version of the nonlo- 
cal NJL model proposed in y^. The function g{z) is 
replaced by a coupling constant G times a normalized 
distribution C{z): 



g{z):=^C{z), jd'zC{z) = \. 



(86) 



The standard (local) gauged NJL model follows for the 
limiting case C{z) = 8'^{z) with J d'^zC{z) — 1. 

In contrast to |22l |. however, G and C are determined 
in conjunction with the behavior of the Polyakov loop L 
or if at temperature T: using the Fourier transform g(j)) 
of g, they are expressed as 



G 

2" 



gip = o)., c{p)^gip)/g{p = o), (87) 



where 



Gip) 



g^ F^{p) + F^^ip) + Fo{p) 



p2 + (T^)2 + 2Tippo 



Note that Fip{p = 0) and hence G diverge at T = 0. This 
comes from an improper treatment of the T = part. 
To avoid this IR divergence at T = 0, we add the T = 
contribution Mq ~ Af|- and replace Fip{id) by 



F^itd) 



1 



and 



1 

' (za^)2 + (r(p)2 + 2T^tdo + M§ ' 



^^^^^ p^ + {po + T^r + Mi 
1 



(89) 



(90) 



F. 



i'i-,n5T) 



(84) 



In fact, such a contribution |Mq comes in G^^ as an 
additional term M§S^^ from the 0{^^) terms (Note 
that 0{^^) terms are absent for G = SU{2)), as already 
mentioned in the above. 

Another way to avoid this IR divergence is to introduce 
the regulator term which is needed to improve the one- 
loop perturbative result and obtain a nonperturbative 
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one according to the Wilsonian renormalization group: 



d^^P 1 ^. 



(27r)^ 2 



:r^^i-pW''g,.Rkip)]3rf{p), (91) 



where fc is the RG scale and Rk{p) is the Fourier trans- 
form of Rk{id). The regulator function R^ introduces a 
mass proportional to k^ , which plays a similar role to Mq 
in the above, as lon g as fc > 0. 

The NJL model [l^ is well known as a low-energy 
effective theory of QCD to describe the dynamical 
breaking of chiral symmetry in QCD (at least in the 
confinement phase), see e.g. [l^ [lj|. The the- 
ory given above by S*^ = S^g"^ + S^g is able 
to describe chiral-symmetry breaking/restoration and 
quark confinement/deconfinement on an equal footing 
where the pure gluon part 3^^° describes confine- 
ment/deconfinement transition signaled by the Polyakov 
loop average. We can incorporate the information on con- 
finement/deconfinement transition into the quark sector 
through the covariant derivative V^f^] and the nonlocal 
NJL interaction Q {G and C), in sharp contrast to the 
conventional PNJL model where the entanglement be- 
tween chiral-symmetry breaking/restoration and confine- 
ment/deconfinement was incorporated through the co- 
variant derivative VlY] alone and the nonlocal NJL in- 
teraction Q is fixed to the zero-temperature case. In our 
theory, the nonlocal NJL interaction Q {G and C) is au- 
tomatically determined through the information of con- 
finement/deconfinement dictated by the Polyakov loop 
L (non-trivial gluon background), while in the nonlocal 
PNJL model (22| the low-momentum (non-perturbative) 
behavior of C was not controlled by first principles and 
was provided by the instanton model. 

To study chiral dynamics, it is convenient to bosonize 
the gauged nonlocal NJL model as done [22]. The nonlo- 
cal gauged NJL model (I55t can be bosonised as follows. 
Define 



$a(a;) := (cr(a;), 7r(a;)). 



(92) 



To eliminate the quadratic term in the nonlocal currents, 
we insert the unity: 



1 = 



/ X>crX>7rexp| - / d^zCiz) 



X / d^a;— [$„(x) + G4>{x + z/2)TMx - z/2)] 
X [$„(a;) + Gi^ix + zl2)T^il^{x ~ z/2)]*}, (93) 

where we have used J d'^zC{z) = 1. Then we have the 
gauged Yukawa model: 



PV^Ve" "«' 



gYi 



VipViP / VaV7fexp{-S^g} 



(94) 



where with x' :^ x + z/2, y' :— x — z/2, 



S: 



«Y 



g - I d^x' I d*2/V(x') 



d\x' - y'){-ij'^V^[y] +mg + Z7V«) 



+ lc{x'-y')r^[<^^{^^ 

+ j d^x^^^{x)<^l{x), 



kC-^)] Hy') 



(95) 



-'off 



d'^p d'^p' 



■V'(p) 



(27r)4 (27r)4 
{2tt)^S\p - p'K-riP, + 9Mp)) +^, + *7V,) 

+ lci^^)r^[<^^{p-p') + -$:ip~p')]\^ip') 



d*p 1 
(27r)4 2G 



'^aipWaiP)- 



(96) 



Finally, the bosonized theory of the gauged NJL model 
is obtained by way of the gauged Yukawa model by inte- 
grating out quark fields as 



QgNJL 






(97) 



where the bosonised action 6*^0^°" is 



5X™ = -Trln 



{S\x' - y')i-ir{d^ - W^^) + nV,) 
+ m, + \c{x' - y')r.[$„(^^) + *:(^^ 



)]} 



J d^x^<^^{x)<^l{x) 



(98) 



S. 



boson 
eff 



Tr In {(2^)454 (p _ p')[-7^bp + 5^^) + ^Vg] 
m, + ^CiP±^)r^[<^^ip-p') + <^lip-p')]] 



d'^p 1 
(27r)4 2G 



^Up)Kip)- 



(99) 



VIII. IMPLICATIONS OF THE POLYAKOV 

LOOP FOR CHIRAL-SYMMETRY BREAKING 

AT FINITE TEMPERATURE 

The thermodynamics of QCD can be studied based on 
our effective theory derived in this paper in the similar 
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way to the nonlocal PNJL model [221. But this must be 
done by including the effect of giuon properly. In the 
PNJL model, the effect of the gluon was introduced by 
the standard minimal gauge coupling procedure, i.e., re- 
placing the normal derivative d^ by the covariant deriva- 
tive D^[£^] :— dfj, — igsi/fj,. In our effective theory, the ef- 
fect of the gluon is introduced through the NJL coupling 
constant G and the nonlocality function C, in addition 
to the minimal coupling i?^['^]. The nonlocal NJL in- 
teraction among quarks are mediated by gluons at finite 
temperature in QCD. Therefore, both G and C charac- 
terizing nonlocal NJL interaction inevitably have tem- 
perature dependence, which could be different depending 
on whether quarks are in confinement or deconfinement 
phases. 

At T = 0, QCD must be in the hadron phase where the 
chiral symmetry is spontaneously broken, which means 
that the NJL coupling constant G{0) at zero temperature 
must be greater than the critical NJL coupling constant 

Gr-. 



G(0) 



, 1 

M2 



>Gc. 



(100) 



The nonlocality function or the form factor C{p) atT = 
behaves 



^^P^ = Y^TJ^^^^^^ = Y^^ 



and 



C(P) 



M§ 



p^ + M§ ■ 



(101) 



(102) 



As an immediate outcome of our effective theory, this 
determines the temperature-dependence of the coupling 
constant G of nonlocal NJL model. Using ([57]) . we have 



G{T) = ig^ 



(Tipr 



Mi 



1 

Ml 



(103) 



which lead to the NJL coupling constant normalized at 
r = 0: 



G(T)/G(0) 



2Af2 



[T^y 



Mi 



1 



(104) 



In the presence of the dynamical quark niq < oo, the 
Polyakov loop is not an exact order parameter and does 
not show a sharp charge with discontinuous derivatives. 
Even in this case, we can introduce the pseudo critical 
temperature TJ as a temperature achieving the peak of 
the susceptibility. Below the deconfinement temperature 
T*, i.e., T < T*, therefore, L ~ or (^ ~ tt, the NJL 
coupling constant G has the temperature-dependence 



G(r)/G(0) ^ i 



2Mi 



7r2T2 



Mi 



+ 1 



(T < TJ 



(105) 



This naive estimation gives a qualitative understand- 
ing for the existence of chiral phase transition. Since 



G{T) is (monotonically) decreasing as the temperature 
T increases, it becomes smaller than the critical NJL 
coupling constant 



G(0) = g' 



1 

Mi 



>Gc, rtoo=^GiO. 



(106) 



Thus, the chiral transition temperature T^ will be de- 
termined (if the chiral-symmetry restoration and con- 
finement coexist or the chiral symmetry is restored in 
the confinement environment before deconfinement takes 
place, i.e., T^ < TJ ) by solving 



GiW 



G(0) 



2Mi 



n^' 



Mi 



Ge 



(107) 



Here we have assumed that the nonlocality function 
C(p) gives the dominant contribution at p = 0, namely, 
C{p) < C(0) = / d^zC{z) — 1 and that the occurrence of 
the chiral transition is determined by the NJL coupling 
constant alone. 

At finite temperature T, the form factor reads 

C{x-y) 



-^E 



nez 



(27r)3 

_fp_Ml 

(27r)3 ~3~ Lp2 + (a; + T^)2 + Mi 



+ 



Jp-{x-y) 



p2 + a;2 -^ Mi 
<Pp Miy_ sinh(ep/T) 



(27r)3 &tp L cosh(ej,/r) — cos((/?) 
sinh(ep/T) 



Jp-{x-y) 



cosh(ep/T) - 1 
where we have defined Cp := ^/^ + Mi and used 
1 _ 1 sinh(ep/T) 



(108) 



^E 



/_. (^ + c)2 + el 26p cosh(ep/T) - cos(G/T) ' 

(109) 

The form factor C does not change so much around the 
deconfinement temperature T ^ T^ (or tp ^ tt). This is 
reasonable since the form factor is nearly equal to the ^ 
correlator Q^^, as already mentioned in sec. II. 

For more precise treatment, we must obtain the full 
effective potential Voft (c, ^) as a function of two order 
parameters a (or {'ipip)) and (p (or (L)), and look for a 
set of values {<J,ip) = (ctojVo) at which the minimum 
Vcs{o'o,'Po) of Vcsi<^,f) is realized. Then if must be re- 
placed by (po in the above consideration. For this goal, 
we must develop the RG treatment for the full theory. 
This issue will be studied in a subsequent paper. 
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IX. HOW TO UNDERSTAND THE 

ENTANGLEMENT BETWEEN CONFINEMENT 

AND CHIRAL SYMMETRY BREAKING 

To discuss the entanglement between confinement and 
chiral symmetry breaking, we wish to obtain the total 
effective potential yQ'-'^ of QCD written in terms of two 
order parameters, i.e., the Polyakov loop average (L) and 
chiral condensate (V'V') , so that its minima determine the 
vacuum for a given set of parameters rriq , T and fj,q when 
Nc and Nf are fixed. Here iriq f oo is the pure Yang-Mills 
limit and rUq J, is the chiral limit. 

The effective potential for the quark part is obtained by 
integrating out quark degrees of freedom. The simplest 
form is obtained e.g., from the bosonized model as 

T^quaik 



= - Tr In 



{^rd^ 



Ca - gj^a'^ + ^^lql'^^ + —(J 



2 

(110) 

Then the RG-scale k dependent effective potential V'^^"'*'' 
for the quark part must be given as the solution of the 
flow equation. In the same approximation as the above, 
it is written in terms of two order parameters a and (/s: 



V, 



quark 



riGZ 



(2^ 



tr In 



i7° - P]l' + Tn-q + C{p)a 



TipTsj + i^qj + R 



,quark 



2G 



(111) 



where T3 — (T3/2 and i?^"^"^ is the regulator function for 
quarks. In the limit fc | 0, indeed, V'^^"'^'' (|llip reduces 
to yq"ark pjoi The effective potential v^q"^'^ pTTj) de- 
pends on TOq, T and /ig when Nc and Nf are fixed. Due 
to thepo dependence of the "mass" function M{p) which 
is an immediate consequence of the nonlocality of the 
present NJL model, it is difficult to obtain the closed 
analytical form by performing the summation over the 
Matsubara frequencies. 

In our strategy, a full effective potential V^^ j, (ct, ip) of 
QCD is given by summing three parts: 



vTu{<^, V) - Vt\v) + Vr'\a, V) + 



^^QCD 



with the pure gluon part V^ "'^('p) = Vr^k (E 
Vt'iv) =Tr{ln[G^^ +5^^i?fc]}. 



(112) 



=Tr{ln[-5^^a^ + (5^^ - 5'''6"'){T^f 

+ 2e^^3j.^5^ ^ <5^^i?fc]}, (113) 



the quark part (|llll) . 



V, 



quark 



{cj,^)^—a^^Tv\n{i-i^d, 



2G 



Ca — TfT^'y 



->quark 



+*^,7^-Ki?r'n' (114) 



and a non-perturbative part AVJ7 (cr, ip) induced in the 
RG evolution according to a flow equation. We assume 
that the total effective action of QCD obtained after inte- 
grating out the fields other than those relevant to chiral 
symmetry and confinement is the form 

Tfe = j dxi j d'x[^Za[dj%{x)f + ]^Z„[d,a{x)f 



V. 



QCD 



eff,fc 



(ct,(^)|, 



(115) 



and obeys the flow equation: 



a^Ffc =-Tr . 







.r.A + i?, 


-1 




t V 1 



dtRk 



dtR, 



't-tik 



(116) 



If the flow equation was solved, we would have ob- 
tained the effective potential of QCD, V^ j, {a, ip) which 
has the following power-series expansion with respect to 
two variables a and ip in the neighborhood of the tran- 
sition point where a ~ ~ L according to the Landau 
argument (as demonstrated in the pure glue case). 



K 



^f'''{a,p)^VHp) + V'^{a) + V^ia,p), 



VHP) =Co + C,p + ^p' + ^^3 + ^^^ + 0{^% 



2 " 4 

VHa,p>)=Ficr^p + 



(117) 



where V^{(p) denotes a part written in terms of ip alone, 
and V^la) denotes a part written in terms of a alone, 
while V'^{<T, ip) denotes the cross term between a and ip. 

Once dynamical quarks are introduced, the exact cen- 
ter symmetry in pure Yang-Mills theory is no longer in- 
tact. Therefore, the QCD effective potential includes 
the explicitly center-symmetry-breaking term. For G = 
SU(2), the center symmetry ip — )■ —ip is explicitly bro- 
ken as Ci 7^ 0, C3 7^ in V^ip) and Fi 7^ in V^ia, p). 
The existence of the cross term is important to under- 
stand the entanglement between center symmetry and 
chiral symmetry, as pointed out by [19|. In fact, the 
one-loop calculation leads to Ci — 0.97434iV/ > and 
Fi — — O.IOGIOSA^/ < {fiq = case) which appears to 
be a good indication for this purpose and serves as the 
initial condition in solving the flow equation. 

In the paper by Schaefer, Pawlowski and Wambach 
[26| , a sort of back-reaction from quarks has been intro- 
duced to improve the effective potential of the Polyakov 
loop, while the NJL coupling remains local. In contrast, 
this paper introduces a back-reaction from gluons to im- 
prove the NJL interaction, leading to the nonlocal NJL 
coupling. However, this does not mean that two treat- 
ments are considered to be alternative. In the presence 
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of dynamical quarks, the running coupling a is changed 
due to fermionic contributions. In [2g], this effect has 
been taken into account as a modification of the expan- 
sion coefficient in the effective potential of the Polyakov 
loop, resulting in e.g., the Nf flavor-dependent deconfine- 
ment temperature Td{Nf). Remembering that the input 
of our analysis is just a running coupling, a sort of back- 
reaction from quarks considered in [26] is easily included 
into our framework by using the running coupling modi- 
fied by quark contributions. Thus, the treatment in this 
paper is already able to take into account back-reactions 
from quarks and gluons mentioned above. 

This section is a sketch of our strategy of understand- 
ing the entanglement between center symmetry and chi- 
ral symmetry. The detailed analysis will be given in a 
subsequent paper. '^ 



X. CONCLUSION AND DISCUSSION 

In this paper, we have presented a reformulation 
of QCD and suggested a framework for deriving a 
low-energy effective theory of QCD which enables one 
to study the deconfinement/confinement and chiral- 
symmetry restoration/breaking crossover transition si- 
multaneously on an equal footing. A resulting low-energy 
effective theory based on this framework can be regarded 
as a modified (improved) version of the nonlocal PNJL 
model [22|. In our framework, the basic ingredients are 
a reformulation of QCD based on new variables and the 
flow equation of the Wetterich type for the Wilsonian 
renormalization group. 

A lesson we learned in this study is that a perturba- 
tive (one-loop) result can be a good initial condition for 
solving the flow equation of the renormalization group 
to obtain the non-perturbative result. In gluodynamics, 
recently, it has been demonstrated [3l|, |33 that the exis- 
tence of confinement transition, i.e., recovery of the cen- 
ter symmetry signaled by the vanishing Polyakov loop av- 
erage can be shown by approaching the phase transition 
point from the high-temperature dcconfinement phase in 
which the center symmetry is spontaneously broken. In- 
deed, the effective potential for the Polyakov loop ob- 
tained in the one-loop calculation which we call the Weiss 



^ It is known that appearance of a mixed- term a'^ip plays an es- 
sential role in the chiral-confinement entanglement. Such a term 
appears in the original PNJL model and leads to the 2 crossovers 
happening almost simultaneously. In the following paper posted 
to the archive after this paper was submitted for publication, it 
has been shown that an effective Polyakov loop-dependent four- 
quark interaction derived by this paper yields stronger correla- 
tion between the chiral and dcconfinement transitions, making 
T^ ~ Td more tightly, than the usual PNJL model. 
Y. Sakai, T. Sasaki, H. Kouno and M. Yahiro, Entanglement be- 
tween deconfinement transition and chiral symmetry restoration, 
e-Print: a rXiv:1006.3648 [hep-ph]. 



potential leads to the non-vanishing Polyakov loop aver- 
age, i.e., spontaneous breaking of the center symmetry. 

For gluon sector, to understand the existence of con- 
finement transition by approaching from the deconfine- 
ment side, we have given the Landau-Ginzburg de- 
scription in the neighborhood of the (crossover) phase 
transition point by analyzing the flow equation of 
the functional renormalization group. The deconfine- 
ment/confinement phase transition is consistent with the 
second order transition for G = SU{2), while the first or- 
der transition is expected for G = SU{3). The detailed 
study of the SU{3) case will be given in a subsequent 
paper. 

The input for solving the flow equation was just a run- 
ning gauge coupling constant, in sharp contrast to the 
PNJL model including several parameters. From the 
viewpoint of a first-principle derivation, this is superior to 
phenomenological models with many input parameters. 

For quark sector, it is possible to obtain the chiral- 
symmetry breaking/restoration transition from the first 
principle. However, we need more hard works, especially, 
to discuss the QCD phase diagram at finite density and 
the critical endpoint. A possibility in this direction from 
the first principle of QCD was demonstrated in one-flavor 
QCD based on the FRG [77]. It will be possible to treat 
chiral dynamics and confinement on an equal footing 
based on our framework along this line [78[. Still, how- 
ever, we must overcome some technical issues to achieve 
the goal of understanding full phase structures of QCD. 
The detailed studies will be hopefully given in a subse- 
quent paper. 
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Appendix A: Reformulation of QCD 



We apply the decomposition ^ to QCD Lagrangian. 
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The quark part is decomposed according to ([2]) as 

S£q ■.=4){i-1>'V^[s^] - mo + 1^7°) V- 

=i!{i-1^'V^ [f] - mo + iM7°)V' + gi'l'' ^A^^, (Al) 

where the covariant derivative I?^ [f] is defined by 

V,[r]:=d^-igY^. (A2) 

The Yang-Mihs part is treated as follows. For the gen- 
eral decomposition ^(x) = l^(a;) + ^(x), the field 
strength ^^^ is decomposed as 

-ig[^f,,^,], (A3) 

where the covariant derivative Df^[y] in the background 
field Yij is defined by 



D^[r]:=d^l^ig[y„-], 



(A4) 



or, equivalently, 



D^if]^^ := d^S^'^+gf^^'^rf. (A5) 

The Lagrangian density ^ym = -i^^^^K] ' ^'"'K] of 
the Yang-Mills theory is decomposed as 

1 






(A6) 






2 

1 

1 

2 
1 



+ o^p.M•*5[^^^''] 



+ -(D^m^, - i^.m^^) • i5[^^ ^i 



1 



(,*3[=^/17 '^t'Jj 



(A7) 



Here the third term on the right-hand side of the above 
equation is rewritten using integration by parts (or up to 
total derivatives) as 






2 

1 

'2 



{D^[nD,[r])^"}s^ 



-l3r'^^{-{Dp[y]D,[y])^^g,^ + gf^^^^^^[y] 



+ D^[r]^^D4yf^}^'' 



(A8) 



where we have used 



gfABC^C^^y^^ 



(A9) 



Thus we obtain 



^^ 



YM 



_^^^^A^AB^^B 



-{tg[^^,^,])\ 



(AlO) 



where we have defined 



W^f := - {D,[Y]DP[y])''^g,^ + 2gf^^^^^^[Y] 



D, 



]AC 



D,[ 



]CB 



(All) 



In the usual background field method, the O(^) term 
is eliminated by requiring that the back ground field Y 
satisfies the equation of motion D^[y]^'"'[y] = 0: 






) = o. 

(A12) 



In our framework, Y do not necessarily satisfy the equa- 
tion of motion. Nevertheless, the 0{^) term vanishes 
from the defining equations which specify the decompo- 
sition. For G = SU{2), D^,[r]n = and ^^ • n = lead 
to 

=G^''[ap( JT, • n) - ^; • D^[r]n] - 0. (A13) 

In order for the reformulated theory written in terms of 
new variables to be equivalent to the original QCD, we 
must impose the reduction condition j38| : 



This eliminate the last term of W^^ in (|Alip . 
Moreover, the 0{,'^^) term is absent, i.e.. 



(A14) 



\{D^[r].T, - D,{r\a^,,) ■ ig{.9::\ jr^ = o, (ai5) 



since DJ^ 
See [3i,ll|, 



^ - D^{r\SC^, is orthogonal to [JTA', JT" 



51|. 



Thus, the Yang-Mills Lagrangian density reads 



(A16) 
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where we have defined 

(A17) 
For G ^ SU{2), the 0{^^) term is absent, be- 
cause ,^fj,„[y] and —ig[^p,,^v'] are parallel to n (this 
is also the case for the sum ^^,y['^] — ig[^^, ^u]), while 
Df,[y]^^ - D:,[y]^f, is orthogonal to n (which follows 
from the fact n ■ ,^ = 0). For G ~ SU{2), therefore, we 
have 



Then the SU{2) gluon part is rewritten into 



(A18) 



^YM = - liG.A^jf - i jr^^g;^,^ms-'^^ 



1 



M^^,^u]? 



(A19) 



where 

(A20) 

Appendix B: Coefficients in the effective potential 

We expand Vj. f, defined by 



V^ i. ^Vw + 4 



(27r)2 
ln(l-2e"Pcos^ + e-2p)}, 



{ln(l-2e-'=^cos(^ + e-2feT) 

(Bl) 



in power series of (f by using the expansion — cos if 

— cos(7 
follows 



cos(7r + if) = cos((^) = 1 — ^(f"^ + ^(fi'^ + 0{(f^) as 



Vr.k ^Vw + 



In 



'^^ dpp^ 



In 



1 + 2e-'=^ - e-'=^(^2 ^ _e-fcT^4 ^ ^~2k^ _^ ^(^6) 



1 + 2e-P - e-Pip^ + —e-P(p'^ + e^^P + o{(p^) 



=H, 



'^^ dpp^ 



PP [ 

7r2 I 



In 



(1 + e-'=^)2 - 6-'=^^ (^2 _^ _e-fcT^4 _^ 0(^6) 



In 



(1 + e-Pf ~ e-Pf + Y^e-V + 0(^«) 



-"IV 



'^^ #p2 



^{>„,i 



-,-fcT\2 



In 



.-fci 



1 - 



^' 



-kT 



ln(l + e-Pf - In 



r^^ 



(l + e-P)2^ 12(l + e'P)2 



(l + e-fer)2 12(1 + 6-'=^)^ 



^' + 0(^«) 



(B2) 



By using log(l + x) ^ x — -^x"^ + 0{x^), therefore, Vj, f. has the polynomial expansion 






(B3) 



where the coefficient is given by the integral form: 



A^ 



2,k 



2 



^-^ dpp^ 



-kT 



4! 127r2 



'^^ dpf 



(l + e-P)2 {i + e~kT)2 
-6e-2p + e-P{l + e-Pf -6e-2feT _^ £-'=^(1 + e-'=^)2 



^o.fe 



'^^ d#2 



i^{'°<' 



-feT\2 



12(l + e-P)4 
-ln(l + e-P)2}. 



12(1 + 6-*:^ 



feT\4 



(B4) 
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The integration can be performed analytically and the coefficient has the closed form: 



A 



2,fe 



3(l + e^)' 



21og(l + e^),s-2Li2(-e^)-y 



A 



i,k 



1 e^" (-2s3 + (s^ + 6) cosh(s)s + 6s + 3 (s^ - 2) sinh(s) - 3 sinh(2s)) 



4! 



127r2 



18(l + e«)V2 
I 



s—kj 



(B5) 



where Li„(z) ~ PolyLog[n, z] is the polylogarithm func- 
tion, and in particular, the dilogarithm satisfies Li[2{z) ~ 



Iz ^"'^K ^' ^^ which is known as the Spence integral. 

is monotonically de- 



fO log(l-t) ^ 

Note that the function r-^f^-^',2 
creasing in x and hence the second term in A2^k is posi- 
tive (non- negative) . The coefficient A2^k is negative and 
monotonically increasing in k and approaches zero for 

k — > oo. 



A. 



2.fe 



1 

6' 



-g<^<0 for fce[0,oo), (B6) 



1 _d^y 
3 - 9^2 ^Tl 



ip=0 



<0 for fce[0,oo). (B7) 



This is because 



kT 



''''JTT^ 



[1 + e py 6 



kT g-fcj 



dpp 

(1 + e 



" 3 ^(1 + 6-^^)2 



(B8) 
(B9) 



The coefficient ^4^^ is positive and approaches for 
fc — > 00, although A4^^k is not monotonically decreasing 
in k. 



^4.0 



1 



A. 



A.k 



4! 127r2 

This is because 

f*kT 



4! 



>0 for /ce[0,oo). (BIO) 



dpp 



dpp 



-Qe-^P + e-P{l + e-Pf 



12{l + e-PY 



1 
12' 



— 6e 



-2feT 



-,-kj 



(l + e-''^) 



-kT\2 



-kT\i 



-6e 



-2kj 



12(1 +6-'=^) 

■(1 



12(1 



-kT\i 



(Bll) 
0. (B12) 



Appendix C: Flow equation for the coefficient 

Suppose that AF^, is of the form: 

(CI) 



The left-hand side of the flow equation reads 

+ 9^^^^ + 0(^«). (C2) 



'^ 3! 



-LP 



The flow equation df,an,k for the coefficient of cp" is ex- 
tracted by differentiating both sides of the flow equation 
n times and by putting cp — 0. The left-hand side is 



9" 



Q^n k 



0=0 



Define 



/(^) 



Anai 



iVT.k+^Vk)- 



(C3) 



(C4) 



The right-hand sides of the flow equation dj^an^k are cal- 
culated from 



_d_ 

dip 



1 



9(^2 



93 



1 + dlfi^) 
1 



-dlfiv) 



1 + 92/(^) 



[1 + 92/(^)]2 



-5^/(^) 



(C5) 



l + dlli^) 



[1 + 92/(^)]2 

[l+a2/(^)]3' 

"g^/(y) 

ll + 92/(<^)]2 

-39^/(y>)9^/(y) 



(C6) 



6 



[l+a2/(^)]^ 



(C7) 



d(p^ 



1 



l+a2/(^) 



-dlfiv) 



[1 + 92/(^)1 



+ 6- 
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-Adlf{^)dlf{^) - 3[9^/(y)] 
[l + a2/(^)]3 

[l + 52/(^)]4 
[l + e/(^)]^'"' 



(C8) 
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If / is an even polynomial in (p, then the flow equation 
is simplified: 






d_ 

dip 

dip'^ 



1 



i + a^/M 



1 + dlfi^) 



-dUiv) 



dtas 



fe"3,fc 



93 



9(^ 



n3 



1 



5^«4,. ^ ^ 



1 + d^fi^) 
1 



1 + dUiip) 



-d'fi^) 



[1 + 5^/M]^ 



[l + dlfM^ 



ip=0 



ip=0 



0, 



(C9) 



ip=0 



ip=0 



ip=0 



ip=0 



= 0, 



(CIO) 
(Cll) 



ip=0 



(C12) 



Therefore, with an initial condition, ai^k = = a^k at 
k = A, the flow equations in the above 



9fcai,fc = 0, dj,a3,k = 0. 



(C13) 



guarantee the solution 



ai,fe=0, a3,fc = (0<A:<A). 



(CM) 
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